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We discuss Kraan - van Baal - Lee - Lu (KvBLL) solutions of the classical Yang-Mills equations with temperature 
in the context of SU(3) lattice gauge theory. We present discretized lattice versions of KvBLL solutions and other 
dyonic structures, obtained by cooling in order to understand their variety and signature. An analysis of the 
zero modes of the lattice Dirac operator for different fermionic boundary conditions gives clear evidence for a 
KvBLL-like background of finite T lattice subensembles with Q = ±1. Using APE-smearing we are able to study 
the topological charge density q(x) of the configurations and to corroborate this interpretation. 



1. INTRODUCTION 
SPECT TO SU{2) 
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The discovery of new caloron solutions with 
non-trivial holonomy has revived hopes for a 
semiclassical description of Yang-Mills theory at 
T 7^ 0, in particular of the deconfinement transi- 
tion. A new aspect of these solutions is the inter- 
play between topological charge and the local and 
asymptotic Polyakov loop. One new feature, the 
separation of a Q — 1 object (caloron) into N c 
dyonic (D) constituents, happens only in some 
part of the parameter space which seems to be 
statistically enhanced in the confinement phase. 

For a couple of years some of us have been 
searching [2] for this type of semiclassical back- 
ground in SU(2) lattice Monte Carlo configura- 
tions at finite T, mostly using cooling. Besides of 
opening a zoo of possible lattice solutions, a re- 
markable outcome was that, near T c , the confine- 
ment phase is distinguished from deconfinement 
by a finite yield of non-trivial caloron configura- 
tions with 2 dyonic constituents. Sometimes they 
are separated into peaks of action, identifiable 
also by different localizations of the fermionic zero 
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mode for periodic/antiperiodic temporal b.c, but 
always monopole pairs can be seen inside, such 
that the Polyakov loop P(x) changes from +1 to 
— 1 . Somewhat unexpectedly, also DD pairs have 
been obtained by cooling as metastable solutions, 
with action S « Si ns t and equal-sign Polyakov 
loop. This happens only in the confined phase 
where the DD pairs originate from caloron-anti- 
caloron structures. Recently we have repeated 
the SU{2) cooling studies at lower temperature. 
We found now that all calorons consist of overlap- 
ping action lumps, yet the Polyakov loop exhibits 
their non-trivial caloron nature. 

2. SIGNATURE OF SU(3) CALORONS 

Let us now come to SU(3). We parametrize as 
Poo = diag(e 27rVl ,e 27rV2 ,e 27r ^ 3 ) the asymptotic 
holonomy, with [i\ < \i2 < /i3 < M4 = 1 + Ml and 
Ml + M2 + Ms = 0- Let yi, y% and y% be three 3D 
position vectors remote from each other. Then a 
KvBLL pQ caloron consists of 3 lumps carrying 
the instanton action split into fractions (12 — Mi, 
/i3 — fi2 and /Z4 — /X3, concentrated near the j/j. 

We discretized continuum calorons and applied 
a few cooling iterations to adapt it to torus 
boundary condition for the gauge field. Then 
an action S 1.3Si ns t and a topological charge 
Q = ±1 arc found. In Fig. 2] we show two ex- 



Figure 1. Topological charge density in the plane 
through all 3 constituent positions for a nontrivial 
(l.h.s.) and a trivial (r.h.s) KvBLL caloron. 

treme cases of lattice KvBLL solutions with well 
separated positions yc a non-trivial one with al- 
most equidistant phases [ii and a caloron of trivial 
holonomy, i.e., Too » diag(l, 1, 1), /ii,/i 2 ,M3 ~ 0. 
In the trivial case one of the constituents carries 
(almost) all action and topological charge. 

In Ref. |2| the zero mode ipo of the Dirac oper- 
ator in the background of a KvBLL solution was 
analytically computed, introducing an arbitrary 
phase in the temporal b.c. of the Dirac operator, 
i.e., ipo (t + 1 /T, x) — exp(2ni£)ipo{t,&)- An in- 
tricate interplay between the phase parameter £ 
and the holonomy phases fa was found. The zero 
mode is not localized simultaneously at all dyonic 
constituents. Instead, it chooses one of them and 
changes its position with changing value of £. A 
selection rule states that the zero mode is local- 
ized at the position jji of the i-th constituent if 
fa < £ < m+i. This implies that for a set of well- 
separated fa and yi the zero mode will visit all 3 
constituents during one cycle of ( from to 1. 3 

We demonstrate this effect using the chirally 
improved lattice Dirac operator 0] . This operator 
has good chiral properties and is known to repro- 
duce the zero mode of latticized standard instan- 
tons |3j . We clearly reproduce the behavior of the 
zero mode for the non-trivial caloron (see Fig. [2}. 
We plot the scalar density of the zero mode. This 
gauge invariant observable is obtained by sum- 
ming |-0o(^c) | 2 at each lattice point x over color 
and Dirac indices. In order to assess the phys- 

3 For cooled SU(2) configurations evidence for KvBLL- 
type configurations was obtained from periodic and anti- 
periodic zero modes of the Wilson-Dirac operator El. 



Figure 2. Scalar density for the non-trivial 
KvBLL caloron of Fig. 1. We show 2 zero modes 
obtained with Q — 0.5 (l.h.s.) and ( = 0.8 (r.h.s.). 

ical relevance of calorons, one should study the 
correlation between the zero mode positions and 
the topological charge density q(x) and find out 
the statistical distribution in the caloron parame- 
ter space. This could be difficult because generic 
equilibrium configurations might contain dyonic 
constituents difficult to separate. Yet it is indis- 
pensable because the calorons have to be disen- 
tangled from the additional topological structure 
which might be of different nature. 

3. COOLED CALORONS AND OTHER 
OBJECTS 

In order to explore the appearance of more 
general SU(3) KvBLL solutions on the lattice 
we cooled hot configurations in the confinement 
phase on 16 3 x 4 lattices (created with Wilson 
gauge action at /3 = 5.5), keeping torus b.c. 
for the gauge field. Then the evolution of the 
holonomy, starting from an average Polyakov loop 
P « (where P — gtrP), is not predictable. 
With increasing number of cooling steps n coo i the 
asymptotic Polyakov loop (an average over 
all lattice points where the local action density 
s < 10~ 4 ) becomes more and more scattered over 
the complex region accessible for the Polyakov 
loop. One can find cooled configurations with any 
non-trivial holonomy. Quasiclassical lattice con- 
figurations are identified when the cooling history 
S(n coo i) goes through a plateau i.e., the equa- 
tions of motion are minimally violated. The dis- 
tribution in Pqo is more centered around the ori- 
gin for higher plateaus and more spread for the 
caloron plateaus S — Si ns t- Of course, only part 
of the caloron events with asymptotic holonomy 
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Figure 3. Polyakov loop distributions over the 
volume for the constructed non-trivial caloron of 
Fig. 1 (l.h.s) and for a caloron obtained by cool- 
ing (r.h.s.). Constituents are marked with stars. 

^oo ^ have such a clean DDD structure 

of the action density as the non-trivial example 
shown in Fig. Q Most of the S = Si ns t plateaus 
have only two distinguishable lumps of action. 
Nevertheless, all caloron events are characterized 
by a distribution of the Polyakov loop field P(x) 
as shown in Fig. [3] where the ideal non-trivial 
caloron is compared with a generic caloron (two- 
lump) event. The bulk of lattice points represents 
the asymptotic holonomy while the constituents 
are mapped to the envelope of this plot. In the 
metric of P, the distance of a constituent from 
the bulk is proportional to the fraction of action 
(topological charge) it carries. 

Purely (anti-)selfdual configurations are most 
frequently observed on the S = Si ns t and S = 
2Si ns t plateaus. The latter, analogous to the 
calorons, preferentially contain four lumps of ac- 
tion, but probably belong to a wider class of Q = 
±2 KvBLL solutions. Besides of these, we have 
found metastable cooled configurations which are 
not of KvBLL caloron type. In the Q — sector 
we found DD events at plateaus S « 0.655j ns < as 
well as DDDD events at S « 1.35Si ns t, usually 
left over from the annihilation of D's and -D's, for- 
mer constituents of caloron and anti-caloron in a 
confinement configuration. 

4. ZERO MODES FOR Q = ±1 EQUILIB- 
RIUM CONFIGURATIONS 

In the last few years it was understood that the 
low lying eigenmodes of the lattice Dirac opera- 
tor provide a good filter which removes short dis- 



tance fluctuations of the gauge field allowing the 
study of topological objects on the lattice. This 
mostly concerns the near zero modes. The pe- 
culiar behavior of the zero mode known for the 
KvBLL configurations suggests that this £ de- 
pendence for any zero mode might be an excel- 
lent, simple and stable tool to highlight the pres- 
ence of caloron background fields in equilibrium 
configurations among other carriers of topologi- 
cal charge, without cooling. A more detailed ac- 
count of the results presented in this section can 
be found in Ref. [S]. 

We have computed eigenmodes of the chi- 
rally improved lattice Dirac operator 0], using 
quenched SU(3) gauge ensembles generated with 
the Liischer-Weisz action [7] on 20 3 x 6 lattices at 
inverse gauge couplings of j3 = 8.20 and (3 = 8.45. 
These values correspond to lattice spacings of 
a = 0.115 fm and a = 0.094 fm and describe 
physics just below and above the QCD phase 
transition (located at j3 c = 8.24 [S]). From the 
full ensembles of gauge configurations we have 
selected those configurations which have only a 
single zero mode, i.e., configurations with topo- 
logical charge Q = ±1. This choice leads one to 
consider only the simplest case where no mixing 
between different zero modes can occur. All or 
a certain fraction of these configurations might 
contain the excess topological charge ±1 in the 
form of well-separated dyons. 

For an equilibrium configuration, let us illus- 
trate the effect of changing fermionic b.c. on the 
scalar density of the zero mode i/'o- In Fig. 01 
we show x-y slices of the scalar density of the 
zero mode for the configuration No. 125 from the 
confined ensemble and compare the cases of anti- 
periodic and periodic temporal fermionic bound- 
ary conditions. The position of the zero mode 
is seen to change under a switch from ( = to 
( = 0.5 which gives us 2 out of 3 expected loca- 
tions. Later we will show that in all 3 positions 
the zero mode is located on top of lumps of topo- 
logical charge with appropriate sign. In the par- 
ticular example shown in Fig. 01 the two positions 
are 1.3 fm apart. 

For a more systematical study we computed 
the zero mode 10 times on the same configura- 
tion using values of £ = 0.0,0.1,0.2, ...,0.9 in the 



Figure 4. Scalar density of the zero mode of 
config. No. 125 (confined phase) for anti-periodic 
(l.h.s.) and periodic b.c. (r.h.s.). In each case we 
show the x-y slice containing the maximum. 

fermionic b.c. This detailed investigation was 
done for 10 configurations in the confined phase 
(/3 = 8.20) and 10 configurations in the decon- 
fined phase [fi = 8.45). 

In the confined phase the expectation value of 
the Polyakov loop vanishes, suggesting a choice 
Mi = -1/3, = 0, n 3 = 1/3, /i 4 = 2/3. 
The selection rule then implies that £ can be 
contained in three different intervals [0,1/3], 
[1/3,2/3], [2/3,1], each selecting another dyonic 
constituent. For 5 of our 10 configurations we 
found that the zero mode indeed visits exactly 3 
different lumps. For 3 of the 10 configurations 
we found only 2 distinct lumps. The explanation 
might be that two of the constituents are very 
close to each other in space. Finally, for 2 of the 
10 configurations we found even 4 different posi- 
tions of the zero mode. This latter observation 
could be an effect of large quantum fluctuations 
on top of the infrared structure mimicking an ad- 
ditional peak. It would be interesting to repeat 
the calculation on a finer lattice with the same 
temperature where one expects the disturbance 
from quantum fluctuations to be less severe. 

In the deconfined phase the Polyakov loop has 
a non-vanishing expectation value. Due to the 
center symmetry of the gauge action this expec- 
tation value can be realized with three different 
Z(3) phases 0, ±27r/3. For real Polyakov loop one 
has hi w (j,2 ~ /X3 f» and ^4 « 1. Roughly 
speaking, the boundary phase parameter £ can 
be contained only in the two intervals [0, 0] and 
[0, 1]. Only for ( w the zero mode is not able 
to see the same lump that it sees for all other 



Figure 5. Histograms of the distance between 
the peak positions of the periodic and the anti- 
periodic zero mode in the confined (l.h.s) and de- 
confined phase (r.h.s.). 

C 7^ 0, which contains a well- localized topological 
charge (Q « ±1). For configurations with com- 
plex Polyakov loop (phase ±27r/3) this critical 
value of £ is shifted to 1/3 respectively 2/3. Our 
set of 10 configurations in the deconfined phase 
contained configurations with real Polyakov loop 
as well as with complex Polyakov loop (here we 
also computed the values £ = 1/3,2/3). For all 
10 configurations we confirmed that only for the 
critical value of £ a different position of the zero 
mode may be taken. 

For a larger set of 89 deconfined configurations 
and 70 confined ones we followed the numerically 
cheaper approach of comparing only the periodic 
(C = 0.0) with the the anti-periodic (£ = 0.5) 
zero mode. For both boundary conditions we de- 
termined the position of the corresponding peak 
and then computed the Euclidean distance 04 be- 
tween the two peaks. The histograms of distances 
are shown in Fig. In the confinement phase 
(l.h.s.) d,4 is distributed over all distances. The 
histogram for deconfinement (r.h.s.) shows that 
for many configurations the mode does not jump 
at all or only over a small d±. There is only a 
small, separate component in the histogram with 
d.4 > 9a, all related to configurations having a real 
Polyakov loop P. We highlight this subsample 
of the deconfined phase in the r.h.s. plot. This 
observation confirms the selection rule for zero 
modes in a KvBLL background. It allows jumps 
between periodic and anti-periodic b.c. only in 
the real sector. This does, however, not imply 
that the zero mode in the real sector always jumps 
over large distances since the two constituents 
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might be close to each other. In Ref. [5| sev- 
eral other tests of the zero modes were performed 
such as a study of their inverse participation ra- 
tio, a measure of their localization, as function 
of C- Throughout we found fair agreement with 
the properties of zero modes in the background 
of a KvBLL solution in the confined phase and 
even excellent agreement in the case of the de- 
confined phase. We have also extended our study 
of zero modes with general boundary conditions 
to lattices with zero temperature geometry and 
still find jumping of the zero modes. A detailed 
writeup of these observations is in preparation. 

In the meantime KvBLL solutions have been 
generalized to higher topological charge, and their 
corresponding zero modes have been constructed 
analytically jHj. It would be interesting to per- 
form a study of the zero modes in higher topologi- 
cal sectors. We also remark that in Ref. ^Oj it was 
argued that the effect of jumping zero modes can 
hardly be due to overlapping standard instantons 
and therefore this effect is truly a consequence 
of the presence of a KvBLL background in the 
Q = ±1 configurations under study. 

5. COMPARISON WITH SMOOTHING 

We have experimented with cooling techniques 
in order to study the correlation between the 
topological charge density and the positions occu- 
pied by the zero mode with changing £. We tried 
restricted cooling [TI], adapted to the Liischer- 
Weisz action, and APE smearing. The latter pro- 
cedure was used by A. Hasenfratz et al. ^2 to 
disclose the topological content of lattice fields. 
In the confinement phase the emerging structures 
are similar for both procedures (which therefore 
can be tuned to each other) and relatively stable 
with respect to cooling/smearing iterations. 

Here we concentrate on results obtained by 
APE smearing with N smear — 10 iterations and 
a = 0.45 As an example for the confined 

phase we show configuration No. 125 on the l.h.s. 
of Fig. [(J The map of the topological density is 
shown in three y-z slices, selected such that each 
contains one of the 3 zero mode positions. Crosses 
show the projection of all 3 positions on the given 
plane, while crosses with circles emphasize the 
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Figure 6. Topological charge density in y-z 
planes, each through one of the zero mode po- 
sitions (marked by the circle): confined config. 
No. 125 with 3 different localizations (l.h.s.) and 
deconfined config. No. 383 with 2 different local- 
izations (r.h.s.). 

zero mode sitting in the plane. There is a com- 
plex pattern of positive and negative charge den- 
sity 4 while the total charge is Q = — 1. The zero 
mode always settles at same-sign lumps which 
are all well-localized. The topological charge is 
not exactly balanced between the lumps, because 
the largest £ interval corresponds to t zm = 5, 
Xzm, = 14, y zm = 16, z zm = 18 (the position 
in the upper plot). 

As an example for the deconfined phase we 
show configuration No. 383 on the r.h.s. of Fig.|SJ 
The configuration has Q = 1 and a real-valued 
asymptotic Polyakov loop. For almost all £ the 
zero mode is localized on top of the well-localized 
constituent seen in the upper plot. For («0 how- 

4 The line separates regions of positive and negative q(x). 
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Figure 7. Cumulative topological charge of same 
sign and opposite sign as function of radius R 
around the center of a zero mode (in the confined 
config. No. 125). 

ever, the zero mode appears much less-localized 
in another region of topological density q(x) > 0. 
The topological charge is much less concentrated 
there, too, consisting of two lumps with the zero 
mode center in between. 

Obviously, also in the deconfined phase one 
finds a complex topological structure which is 
more sensitive to cooling than in the confinement 
phase. The gap separating the near zero modes 
strongly depends on it. We have observed that 
the gap rapidly closes if we cool too much. 

For the confined configuration No. 125 and for 
the three possible localizations of the zero mode 
we calculated the cumulative topological charges 
of same and opposite sign (relative to Q) as a 
function of the radius of a 4-sphere. The curves 
shown in Fig. [7] are an average over the three 
cases. The radius inside of which the topological 
density has a coherent sign of q{x) is a good defi- 
nition of the constituent radius. Then the average 
topological charge of a constituent is Q const ~ ±5 
in the confinement phase. 

6. CONCLUSIONS 

We obtained new classical solutions by cool- 
ing of confined SU(3) lattice configurations. For 
Monte Carlo configurations in the topological sec- 
tor Q = ±1, the dependence of the localization 
of the fermionic zero mode on the aperiodicity 
angle £ matches the specific expectations derived 



from the KvBLL caloron picture for the confined 
and deconfined phases. For slightly APE-smeared 
configurations we are able to corroborate these 
findings by considering the topological charge dis- 
tribution close to the zero mode locations. 
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